Are galaxies extending? 
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It is suggested that the recently observed size evolution of very massive compact galaxies in the 
early universe can be explained, if dark matter is in Bose Einstein condensate. In this model the 
size of the dark matter halos and galaxies depends on the correlation length of dark matter and, 
hence, on the the expansion of the universe. This theory predicts that the size of the galaxies 
increases as the Hubble radius of the universe even without merging, which agrees well with the 
recent observational data. 

PACS numbers: 98.62.Gq, 95.35.+d, 98.80.Cq 

Dark matter (DM) and dark energy [l[ are two of most important unsolved puzzles in modern physics and cosmology. 
Identification of one DM particle species by a direct detection experiment such as LHC or DAMA [2] is not enough 
to fully solve the dark matter problem, because there can be multiple species of DM and we have to explain the 
observed distribution of DM in the universe. Although the cold dark matter (CDM) with the cosmological constant 
(i.e., ACDM) model is popular and remarkably successful in explaining the large scale structure of the universe, it 
seems to encounter problems on the scale of galactic or sub-galactic structures Q . Numerical simulations with ACDM 



O ' model usuallypredict a cusped central halo DM density and too many satellite galaxies compared to astronomical 



observations [J [ESQ 



On the other hand, it is known that the cold dark matter model based on Bose Einstein condensate (BEC) 
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or scalar field dark matter (SFDM) can alleviate these problems pS [HL OS 03 and well explain the observed 
rotation curves of galaxies 0, EH- Henceforth, I will designate the two models as the BEC/SFDM model [HI]. 
In this model the galactic halos are like gigantic atoms where cold boson DM particles are condensated in a single 



macroscopic wave function ip(r). ( The idea of giant atoms as hypothetical stars goes back to Kaup and Ruffini 17 1 
and developed by Schunck and others ( See 1 1811 fo r a review. V) Similar halo DM ideas were sugg ested by many 
authors [13, EL OS HE EL HI US HI HE HUH HI, HI HE HI HI HE Hi, HE HI . (See [H, H3 and references 
therein.) This BEC/SFDM model, now often known as the fuzzy DM model, is a variant of the CDM model. It is 
more about the state of DM particles rather than the DM particle itself. 

The recent observations 13a, l39l lib] ] of the size evolution of massive galaxies even deepen mysteries of DM and 
55 formation of galaxies. In _40], using the combined capabilities of earth-bound telescopes and the Hubble space 
f^*) . telescope, the size evolution of 831 very massive galaxies since z ~ 2 is investigated. In [3i| it is observed that massive 
' quiescent galaxies at z ~ 3 have a very small median effective radius r = 0.9 kpc. According to the observations, very 
• i-H , massive (> IO^Mq) galaxies have a factor of about 5 smaller size in the past (z ~ 2) than their counterparts today. 
■ These compact but massive galaxies are puzzling in the context of ACDM model, because the size evolution is usually 
5-j | attributed to a hierarchical merging process of small galaxies to form a larger galaxy. Thus, it is expected that small 
. - - i early galaxies have small mass too. The finding of compact but very massive early galaxies which have disappeared 
calls this interpretation in question, since these galaxies had almost reached the maximum mass limit observed today 
and might not have been experienced significant merging after z ~ 2 [38| . Various mechanisms depending on visible 
matter could change the size-mass relation by a factor of ~ 2 but not a factor of ~ 5 [3t| . 

In this paper, it is suggested that this galaxy size evolution problem can be also solved in the BEC/SFDM model, 
if the correlation length of the DM condensate is time dependent. First, let me briefly review the BEC/SFDM model. 
In 1992, to explain the observed galactic rotation curves, Sin [S El suggested that galactic halos are astronomical 
objects in BEC of ultra light (mass m ~ 10 -24 eV) DM particles such as pseudo Nambu-Goldstone boson (PNGB). 
In this model the halo cold boson DM particles are condensated in a single macroscopic wave function ip( r ) and the 
quantum mechanical uncertainty principle prevents the halos from self-gravitational collapse, while in usual CDM 
models DM particles move independently and incoherently, ^(r) satisfies the non-linear Schrodinger equation with 
the Newtonian gravity; 
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ihdti* = EiJj = -—V 2 t/j + Vifj(r), (1) 
2m 



where E is the energy of each DM particle. The gravitational potential is given by 

V(r)= f dr'^ f dr"47rr" 2 (p vis (r")+p DM (r"))+V , (2) 
Jo r Jo 

where DM density pdm{t) = M \ip(r)\ 2 and p V i S is the visible matter (i.e., stars and gas) density. M is a mass 
parameter and Vq is a constant. According to the model, the condensation of DM particles of which huge Compton 
wavelength A c = 2nh/mc ~ 10 pc, i.e. m ~ 10~ 24 eV, is responsible for the halo formation. 
The author and Koh 

sin generalized Sin's BEC model in the context of quantum field theory and the general 
relativity and suggested that DM can be described as a coherent scalar field (i.e., SFDM). In this model the BEC DM 
halos are giant boson stars (boson halos [H, S3, H|) described by a complex scalar field <f> having a typical action 

S = J ^~9d A x[^ - ^o;,., : „ - 17(0)] (3) 

with a potential 17(0) = 1J ^\(j>\ 2 + ||0| 4 - In this paper, the case with A = will be considered for simplicity. The 
spherical symmetric metric is ds 2 = —e v ^dt 2 + e x ^dr 2 + r 2 dfl 2 , where r is the radial coordinate. 

In the BEC/SFDM theory, due to the uncertainty principle AxAp > H/2, the characteristic length scale of a BEC 
DM halo £ is inversely proportional to the mass of DM particles, i.e., £ ~ Ax ~ h/Ap ~ h/mAv ~ 10~ 3 c H/m, 
where Av is the velocity dispersion of DM particles and c is the light velocity. This £ should be comparable with the 
characteristic length scale of the DM condensate. (For DM halos or boson stars, the correlation length or de Broglie 
wavelength of the condensate is more suitable for the length scale than the Compton wavelength A c ~ 1/m 0, IH|.) 
From this one can also obtain the minimum mass for galaxies M c = G ^ m i — 10 7 M©, which is recently observed |46| . 

The temperature-dependent correlation length of DM is about the thermal de Broglie length £(T) ~ 2irh c/kT of 
the DM condensate [45| | , which is an increasing function of the time. This leads to a surprising possibility that the sizes 
of a DM halo and a galaxy embedded in the halo are slowly increasing functions of the time even without merging. 
The DM halo provides a potential well to trap the visible galactic matter such as stars and gas. The extension of the 
halo induces the extension of a visible part of the galaxy. To sec this consider the Newtonian limit of the equations 
of motions from the action in Eq. ([3]), or, Eq. (TTJ), which can be written in the dimensionless form as 



V 2 V = (a 2 + p vis ) 
V 2 ct = 2(V-E)a 



2„_onr ^ ■ (4) 




Here a = V^Ge-' lEt iP/c 2 is a dimensionless form of the wave function tp (or the scalar field 0) and r —* mcr/H, 
t — > mc 2 t/h @ and all other quantities are dimensionless too from now on. Since galaxies are dominated by DM, 
we will ignore p V i S and assume that visible matter passively moves inside the potential well V of the DM halo. The 
equations above have approximate solutions [47] 

~ _7fV> fl , (V(0)-E)r 2 \ 

(5) 

To see the size evolution for the most massive galaxies, we need to calculate the DM distribution in their halos for 
a fixed galaxy mass M = J Q 4Trr 2 (a(r) 2 + p V i S (r))dr ~ 0(|cr| 2 r 3 ). Thus, for a constant M and under the scaling 
£ — > corresponding to the increase of the correlation length, the other parameters for DM halos scale as r — > Ir, 
a — > l~^l 2 o and V <~ M/r — * V/l. This means, under this scaling, the gravitational potential well becomes shallower 
and the size of DM halos r^M increases like I. We can also assume that the total energy (kinetic + potential) of a 
star E s t measured from V(0) is conserved during the extension. The star orbits around the galactic center within the 
gravitational potential of the halo V. Note that E st is different from the energy of a DM particle, E. The visible 
radius of a galaxy can be defined by an orbital radius of outermost stars r*, which can be defined as the position 
satisfying the condition E st — V(r = r*), i.e., the position where the kinetic energy of the stars is zero. From Eq. 
and the energy conservation one can see that E st = V(r*) = a(0) 2 r 2 /6 and r* = y/6E st /cr(0). Since cr(0) — > Z~ 3 / 2 cr(0), 
the orbital radius of the star scales as r* — > ^ 3 / 2 r*. It means that the visible radius of galaxy follows approximately 
the 3/2 power of the size of its dark matter halo. (The extension of the halo also reduces the density of matter 
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r(z)/r(0) 
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FIG. 1: (Color online) Observed size evolution r(z) of the massive galaxies versus the redshift z. The red squares denotes 
the size evolution of spheroid-like galaxies and the blue circles are for disc-like galaxies (Data from Fig. 9. of Ref. [13] ). 
The black dot represents the typical compact galaxy at z = 2.3 in Ref. [39f. The line represents our theoretical prediction 
tdm(z) i 'r dm(0) = 1/(1 + z) for DM halos. The dashed line represents the predicted size evolution of the visible parts of the 
galaxies r*(2)/r*(0) = 1/(1 + z) 3//2 which agrees well with the observational data for the disc-like galaxies. 



and the rotation velocity of stars. This is also observed [40(. Since the observed 'extension speed' of the galaxy is 
very low (O(kpc/10 10 yrs)) we can treat the expansion as an adiabatic one. Thus, the temperature is inversely 
proportional to the scale factor R of the universe, i.e., T ~ l/R, as usual.) 

Collecting all together, we obtain a simple relation between the size parameter Tdm{z) of DM halos and the redshift 

z; 

tdm(z) = £(z) = T(0) = R(z) = 1 

r DM (0) C(0) T(z) R(0) l + z U 

Here tdm( z ) denotes the effective radius of a DM halo at z, and rfljvf(O) denotes that at present. Thus, the size of 
visible galaxies r* evolves as 

r_M f r DM (z) \ 3/2 (R(z)\* ( 1 \ f 

r.(0) \r DM (0)J \R(0)J \l + z) ' [ ' 

Since the Hubble radius H^ 1 ^) ~ t and R(t) ~ t% during the matter dominated era, the size of the most massive 
galaxies increases as the Hubble parameter at the first order. Fig. 1 shows the observed size evolution r(z) of the 
massive galaxies [1^, H3| versus the redshift z. The theoretical prediction (r*(z) in Eq. (J7J)) is well coincident with 
the observational data. The black dot represents the typical parameter for 9 galaxies in Ref. [3^| for which I used 
the typical values z = 2.3, r(z) — 0.9 kpc, and r(0) = 5 kpc from the paper. For the error bar for the dot, I used 
the effective radiuses of the smallest and the largest galaxies in the table 1 of Ref. [3^|. Although there are still 
many theoretical and observational uncertainties, the coincidence between the theoretical prediction and the data 
for disk-like galaxies is remarkable, especially for high z. The small discrepancy can be attributed to the systematic 
observational uncertainties and ignoring of possible merging history. It is unclear why disc-like galaxies and spheroid- 
like galaxies seem to show rather different size evolutions. We need more observations to determine whether the 
difference is real or not. 

To be more illustrative, we perform a numerical study using the shooting method 0]. Fig. 2 shows the result of 
our numerical study with boundary conditions dV/dr(0) = 0, V(oo) = and da/dr(0) = 0. We consider 3 cases 
with the parameters (cr(0) = 5 x 1(T 7 ,V(0) = -3.678 x Kr 7 ,£ = -1.52 x 1CT 8 ), (cr(0) = 3.21 x 1CT 7 , V(0) = 
-2.72 x 1CT 7 ,£; = -5.71 x 1(T 8 ), and (cr(0) = 2.17 x 1(T 7 , V(0) = -2.04 x 10~ 7 ,£; = -6.23 x 10~ 8 ), respectively, 
from the top to the bottom for a(0). With E st = 10~ 7 we changed the length scale as I = 1, 1.5,2 for 3 cases and 
obtained r, = 1736, 2861, 4590,respectively. The masses M = 0.013,0.0125,0.0124 are similar for all 3 cases. The 
numerical results support the theoretical argument above. 

We need to check the reliability of the Newtonian approximation used to derive the Eqs. ^ and ([7]). From the action 
in Eq. © and by defining a = y/4irGe~ luJt <f) one can obtain the dimensionless versions of the scalar field equation and 
the Einstein equation [48| , which can be reduced to the equation in Eq. (fT]) in the Newtonian limit @ . Since the typical 
compact galaxies observed have mass M ~ 10 u Mq and size R ~ lkpc, the dimensionless gravitational potential 
V ~ M/R ~ 4.78 x 10~ 6 and the expected rotation velocity dispersion is v rot — 0(W) ~ 0.0022c ~ 650fcm/s, which 



4 



is comparable with the recent observational data v ro t — 510lgQ 5 fcm/s by Dokkum et al (49|. Thus, the compact early 
galaxies are basically non-relativistic objects and the Newtonian approximation is good for these galaxies. 

Schunck, et al pointed out that the effect of pressure generated by scalar fields should be included in the rotation 
curves [l9|, The rotation velocity given by the circular geodesies is v ro t — rdv / dre v /2 ~ M(r)/r + p r r 2 e x+v /2, 
where the second term denotes the contribution of the radial pressure p r = p — U = \{io 2 G 2 er v + o' 2 e~ x — U) of the 
DM field. Since a ~ O(10~ 7 ) <C 1 and v rot «cwe can use the weak field approximation in Ref [5l|, [S^]. A relevant 
parameter for this approximation is e = (1 — w 2 /m 2 ) 1/2 <C 1. In this limit li ~ m, d/dr ~ e and e X ' v — > 1. Thus, 
we obtain (dimensionless) p r ~ \{& 2 + e 2 a 2 — a 2 ) = 0(e 2 cr 2 ) <C p ~ 0(a 2 ) and neglecting the contribution from p r 
m our model (with U = m 2 2 /2) is a good approximation for these galaxies. This is different from the case of the 
model with massless scalar DM particles, where U = fl9l. [Hoj. 
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FIG. 2: (Color online) The dark matter field a (blue dotted lines), the energy level of visible matter (green thick lines) and 
the gravitational potential V (red dashed lines) for a galaxy as a function of distance r from the halo center. As the universe 
expands, the temperature of DM decreases and its correlation length £ increases. This induces the increase of the size of the 
visible galaxy (r*) represented by the green lines. 



We have shown analytically and numerically that the BEC/SFDM theory could explain the observed size evolution 
of the most massive galaxies. In our theory the size of galaxies can increase not only by merging or accretion of small 
galaxies but also by the increase of the length scale of DM halos itself. For small and medium sized galaxies, it is 
hard to distinguish these two effects. This may explain why the pure size extension without merging is observed only 
recently and only for the most massive galaxies. Our theory explains how these massive galaxies were so dense in the 
past and reached the size of the massive galaxies today. 

A self- gravitational redshift effect of the DM halo may contribute to the relation in Eq. © or Eq. |(7J). The observed 
redshift of the galaxies could be a combination of the cosmological (1 + z) and the gravitational redshift (1 + z g ), i.e., 
(l+z)(l+z g ) The gravitational redshift parameter for the boson halo is given by [l9[ z g = eM 00 ) - "^/ 2 — 1 ~ V 
in the weak gravitation limit. Since V ~ O(10 -6 ), the correction from z g is negligible for galactic DM halos in this 
model [HI compared to the cosmological redshift. A time- varying gravitational constant or the gravitational memory 
in Brans- Dicke models could mimic the extension of the galaxies [5^, HH . 

In conclusion, the idea that DM is in BEC seems to provide us a new way to explain not only the CDM problems 
but also the galaxy evolution. From this perspective it is important to determine the exact size evolution of the 
galaxies by future observation. 
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